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ABSTRACT 

We report results from the second and third seasons of observation with the QUaD experiment. Angular 
power spectra of the Cosmic Microwave Background are derived for both temperature and polarization at both 
100 GHz and 150 GHz, and as cross frequency spectra. All spectra are subjected to an extensive set of jackknife 
tests to probe for possible systematic contamination. For the implemented data cuts and processing technique 
such contamination is undetectable. We analyze the difference map formed between the 100 and 150 GHz 
bands and find no evidence of foreground contamination in polarization. The spectra are then combined to 
form a single set of results which are shown to be consistent with the prevailing LCDM model. The sensitivity 
of the polarization results is considerably better than that of any previous experiment — for the first time 
multiple acoustic peaks are detected in the E-mode power spectrum at high significance. 
Subject headings: (cosmology:) cosmic microwave background, cosmology: observations, polarization 



1. INTRODUCTION 

The anisotropy of the Cosmic Microwave Background 
(CMB) gives us direct insight into the structure of the Uni- 
verse when it was a tiny fraction of its current age, and 
is one of the central pillars of the enormously successful 
standard cosmological model. The temperature anisotropy 
power spectrum has now been measured to good precision 
from the larges t angular scales dow n to a small fraction of 
a degree (e.g. Reich ardt et alj d2008l) ) — the expected series 
of acoustic peaks is present and fitting the spectrum yields 
tight constr aints on the basic par ameters of the cosmological 
model (e.g. lDunklev et alj d2008l) ). 

The CMB is expected to be polarized at the ~ 10% level 
principally because of motions in the material at the time 
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of last scattering. Since the plasma flows along gradi- 
ents in the density field the resulting observable polariza- 
tion pattern has gr a dients (E-modes), but zero curl (B-modes) 
(e.g. iHu & White! (Q997)). Given a standard cosmological 
model fit to the temperature spectrum (TT), the E-mode spec- 
trum (EE), and temperature-£'-mode cross spectrum (TE), 
are nearly deterministically predicted — only at the largest an- 
gular scales is there additional information. It is important to 
remember that, although very successful, the standard cosmo- 
logical model (which we will refer to throughout as LCDM), 
contains several components which we have only circumstan- 
tial evidence for (dark matter and dark energy). Measuring 
the EE and TE spectra is thus a crucial test of the overall 
theoretical paradigm. 

As the CMB travels to us through the developing large 
scale struct ure s ubtle deflections due to gravitational lensing 
occur (e.g. IHu! d2003l) ). This converts some fraction of the 
E-mode pattern into the so-called lensing B-modes — this 
effect is most important at smaller angular scales. In ad- 
dition, if the cosmogenic theory known as inflation is cor- 
rect, there must also be large angular scale fi-modes caused 
by g ravity waves propagating through the primordial plasma 
(e.g. lSeliak & Zaldarriagal d 1997b ). 

The polar ization of the CMB was first detected by the DASI 
experiment dKovac et alj [2002). and since then several exper- 
ime nts have reported measurements of the EE and TE spec- 
tra (Barkat s~etaT]|2005l: iReadhead et~al|l2004t [Montrov et alJ 
|2006t ISievers et alj 120071: iPage et alj l2007t , iBiTchoff et alj 
l2008t iNolta et al.ll2.Q08h . Thus far all reported measurements 
are consistent with LCDM, although precision remains lim- 
ited. B-mode polarization has not yet been detected — all 
results so far are upper limits . We previously r eported prelim- 
inary results from QUaD in lAde et al.1 d2008l) . In this paper 
we report considerably improved results using data from the 
second and third seasons of observations. 

This paper is structured as follows: in Section|2]we briefly 
review the instrument and detail the observations, Section [3] 
describes the low level data processing and calibration, Sec- 
tion |4] outlines the steps used to make timestream data into 
maps, Section [5] describes the simulation methodology, Sec- 
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FIG. 1 . — Left: the relative gains of the first ten 150 GHz detector channels as measured by elevation nods every half hour over two seasons. Right: a histogram 
of the percentage fluctuations of the timeseries at left including all channels. 



tion [6] converts the maps into power spectra, Section [7] gives 
the results of jackknife tests, Section [8] describes foreground 
studies, Section [9] gives the final combined power spectrum 
results, Section[lO]contains some further investigations of sys- 
tematic effects, and SectionfTTIstates our conclusions. 

2. INSTRUMENT AND OBSERVATIONS 

The design, implementation, and performance of the 
QU aD experiment is des cribed in detail in a companion pa- 
per dHinderks et alj |2008) hereafter referred to as the "Instru- 
ment Paper" — only a very brief summary will be given here. 
QUaD was a 2.6 m Cassegrain radio telescope o n the mount 
origin ally constructed for the DASI experiment dLeitch et alj 
|2002|) . This is an az/el mount with a third axis allowing 
the entire optics and receiver to be rotated around the line 
of sight (referred to as "deck" rotation). The QUaD re- 
ceiver consi sted of 3 1 p a irs of polarization sensitive bolome- 
ters (PSBs; Uones et all (J2003l) ). 12 at 100 GHz, and 19 at 
150 GHz. The detector pairs were arranged in two orienta- 
tion angle groups separated by 45°. The mount is enclosed 
in a large bowl shaped reflective ground shield on top of a 
tower approximately 1 km from the geographic South Pole. 
The bolometers were read out using AC bias electronics, and 
digitized by a 100 Hz, 16 bit ADC. The raw data were staged 
on disk at Pole and transferred out daily via satellite. QUaD 
was decommissioned in late 2007. 

The observations reported on in this paper were made dur- 
ing the Austral winter seasons of 2006 and 2007 — the QUaD 
telescope was not able to observe during the summer due 
to contamination from the Sun. Complete CMB observation 
runs occurred on 171 days of 2006 and 118 days of 2007 (de- 
fined as an uninterrupted run with the Sun below the horizon). 
Of these available 289 days 44 were rejected after initial low 
level processing — mostly due to very bad weather, with a 
few due to instrumental problems. A further 43 days show ob- 
vious signs of contamination by the Moon, and 59 more fail 
a very conservative Moon proximity cut — Moon contami- 
nation is discussed in the Instrument Paper and Section 110.71 
below. This leaves us with a total of 143 days of data which 
are used in the current analysis. For simplicity the cut granu- 
larity is very coarse — we only consider complete days, and 
if there is anything wrong with any part of a day we cut the 
entire day. It would certainly be possible to include somewhat 
more data with additional work. In Figure Q] we can see the 
resulting set of days used — the monthly gaps are due to the 
rising and setting of the Moon. 

The QUaD telescope is mounted on a tower at one end of 
the MAPO observatory building. At Pole the celestial sphere 
rotates about the zenith every 24 hours. Therefore to minimize 



the potential for contamination from the building, or the heat 
plume from the furnace it contains, each day of observation 
starts at a fixed LST such that our chosen CMB field (centered 
on RA 5.5h, Dec -50°) has just cleared the laboratory build- 
ing. The observations are split into two blocks, each of eight 
hours, with special calibration observations before and after 
each block. Between the two blocks the entire telescope is ro- 
tated by 60° around the line of sight, and then approximately 
30 minutes are allowed for thermal stabilization. The total ob- 
servation schedule takes about 19 hours, with the remainder 
of the 24 hour period being taken up by fridge cycling. 

Each 8 hour block of CMB observations is divided into 
16 half hour periods. Each starts with an observation of the 
internal calibration source followed by an "elevation nod" — 
the telescope is moved up and then down again by one degree 
in elevation injecting a large signal into the data stream due 
to the atmospheric gradient (see the Instrument Paper for de- 
tails). The telescope is then scanned back and forth five times 
over a 7.5° throw in azimuth, with the scan being applied as a 
modulation on top of sidereal tracking. The scan rate is 0.25° 
per second in azimuth translating to around 0.16° per second 
on the sky at our observing elevation. The pointing position is 
then stepped by 0.02° in declination and the process repeated 
four times. Including time for moves and settling, these cal- 
ibration observations, plus the 40 "half-scans" of 30 seconds 
each, take half an hour. 

During the first half hour of each pair the pointing center is 
RA 5.25h. This is then set to RA 5.75h and the declination 
offsets repeated during the second half hour. The observa- 
tions are thus taken in a lead-trail manner where the scan- 
ning pattern is identical in ground fixed azimuth-elevation co- 
ordinates between the two half hours of each pair. By sub- 
tracting the lead and trail data one therefore cancels any sig- 
nal coming from the ground which is constant over the half 
hour, while producing a difference map of the sky. We do 
see significant ground pickup and lead-trail field differencing 
is used thro ughout the analysis presented in this paper (see 
Section[33]). 

From hour to hour the declination offsets are cumulative, 
but they are reset before the second eight hour block. The tele- 
scope therefore scans a "letter box" region 0.64° high in Dec 
twice per day, at two different line of sight rotation angles. 
Each day the observation region steps by 0.64° in Dec to cover 
the entire field as rapidly as possibly (in « 9 days), and then 
cycles around with a 0. 16° offset to generate an eventual four 
fold interleaved coverage pattern (after 37 days). This ob- 
servation pattern was repeated throughout the 2006 and 2007 
seasons a total of 8.5 times. Figure [2] shows the location of 
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FIG. 2. — Location of the QUaD lead-trail fields delineated in white 
on an equal area azimuthal projection about the SCR The color map is 
th e prediction of dust e mission intensity at 150 GHz from FDS model 
8 i Finkbeiner et al. 1999). The color scale is linear from to 100 fiK, and 
is heavily saturated. (The B03 deep and shallow regions are delineated in 
red, while the white asterisk and cross show the locations of RCW38 and 
PKS0537-441 respectively.) 

Note that only azimuth scanning is used. Since the tele- 
scope is only ~ 1 km from the Earth's rotational axis essen- 
tially zero "cross-linking" of the map occurs. For the multi- 
pole range presented in this paper this has only a small nega- 
tive impact on the final CMB power spectrum results. 

The bright galactic HII region RCW38 was observed on 1 1 
days distributed through 2006 and 2007, to monitor the beam 
offset angles and shapes, and the bright quasar PKS0537- 
441 was observed during the 2007 season to further study the 
beam shapes. These observations are used in the analysis be- 
low (see Sections [4.21 and [5.1.1b - 

We also observed several other discrete sources (Cen A, 
Galactic Center, Moon), as well as conducted a survey of part 
of the galactic plane. These observations will be described in 
future papers. 

3. LOW LEVEL DATA PROCESSING 

This section describes the low level data processing which 
occurs before the timestream is binned into maps (Section HJ 
and analyzed to generate simulations (Section[5j. The data is 
deconvolved to remove the effects of detector timeconstants, 
relative gain calibrated, and field differenced. 

3.1. Deconvolution of detector temporal response 

The initial stages of the data reduction are performed on 
each day of data individually. The first step is deconvolution 
of the detector time constants — the bolometers have non- 
instantaneous response to changes in incident optical power, 
and hence the output data timestream is a low pass filtered 
version of the desired input waveform. With knowledge of 
the detector time-constants this filtering can be undone at the 
price of increased high frequency noise. 



We measured the temporal response of our detectors in situ 
using an external Gunn oscillator source as described in the 
Instrument Paper. We find that many of the detectors are well 
fit by a simple single time-constant model. The median pri- 
mary time constant is 16 ms. However a substantial fraction 
(~ 50%) require a second additive time-constant to obtain a 
good fit — we are hypothesizing that some fraction of the in- 
cident heat goes into a second reservoir which has a weaker 
coupling to the thermal bath. In a couple of cases the sec- 
ond time-constant is several seconds long but the two time- 
constant model is still a good fit — we retain these detectors. 
In two other cases the dual time-constant model is not a good 
fit and we reject these detectors leading to the loss of two 
channel pairs (both at 100 GHz). 

To check that the deconvolution process is working, and 
that the timeconstants are stable over time, we examine 
forward-backward jackknife maps of the bright compact 
source RCW38 taken on 1 1 days distributed through 2006 and 
2007. These show no detectable residuals and hence cancel- 
lation to < 1%. 

As part of the deconvolution process the timestream is also 
low pass filtered to < 5 Hz. After deconvolution we de-glitch 
to remove cosmic ray hits etc. (loosing <C 1% of the data). 

For this analysis we also exclude two additional channel 
pairs; one (at 150 GHz) due to a time evolving scan syn- 
chronous signal, and a second (at 100 GHz) due to strongly 
atypical pair differenced noise. We are thus left with 9 of 
the possible 12 pairs at 100 GHz, and 18 of the 19 pairs at 
150 GHz. 

3.2. Relative gain calibration via elevation nods 

We measure the relative gains of the detector channels using 
the elevation nod method mentioned in Section[2] above. The 
airmass through which each channel pair was looking is cal- 
culated from the elevation encoder reading and then regressed 
against the observed signal to yield a calibration factor in volts 
per airmass. We then simply normalize each channel's gain to 
the mean of all channels: 

V/(f) = Vi(f)- (1) 
gi 

where subscript i is a loop over channels, V (t ) are the detector 
timeseries data, and g are the elevation nod gain factors. This 
equalizes the gains both within, and between, channel pairs. 

The nominal accuracy of each elevation nod gain measure- 
ment is -c 1%. Weak trends in the apparent relative gains are 
observed over short and long timescales, the cause of which 
is unknown. However as we see in FigureQ]this leads to fluc- 
tuation over the entire two season time span of only s» 1% 
rms. In this analysis we choose to regard these variations as 
real and the relative gains derived from each elevation nod are 
applied to the subsequent 40 half-scans. 

However even if these apparent variations are false, random 
errors of this magnitude will cause negligible leakage of total 
intensity to polarization (hereafter T to pol. leakage) as, for 
example, T will leak sometimes into +Q and and sometimes 
into — Q, averaging down in the final maps over both time and 
detector pairs. 

Possible systematic errors in the relative gains are of much 
greater concern and are discussed and simulated further in 
Section [103] Imperfect deconvolution would lead to the rel- 
ative gain of a detector pair being a function of temporal fre- 
quency. Since the elevation nods measure the gain at an effec- 
tive frequency well below the CMB measurement band this 
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would result in systematic T to pol. leakage. Using the Gunn 
oscillator derived ti me co nstants and deconvolution procedure 
described in Section [XT] the lack of "monopole" residuals in 
pair difference jackknife maps of the bright source RCW38 
indicate that the low frequency elevation nod gains are accu- 
rate to better than 1% at high (beam scale) frequencies. 

3.3. Examination of timestream and noise spectra 

The sum of the signals from a PSB pair measures the total 
intensity of the incident radiation, while the difference mea- 
sures polarization. Having performed the relative gain cali- 
bration and sum/difference operations Figure [3] shows a typi- 
cal scan set at this stage of the processing. Atmospheric emis- 
sion at 100 and 150 GHz is dominated by oxygen and water 
vapor respectively. The intensity of the radiation received is 
an integral over the temperature and density along the atmo- 
spheric column through which the telescope is looking. Wa- 
ter vapor is poorly mixed in the atmosphere leading to much 
greater fluctuations at 150 GHz. From the lack of an obvi- 
ous scan synchronous component we can infer that the wind 
blows inhomogeneities through the beam faster than the tele- 
scope scans. Since the QUaD beams do not diverge until high 
in the atmosphere the pair sum timestreams are highly corre- 
lated across the array. As expected the atmospheric emission 
is clearly very weakly polarized leading to strong cancellation 
in the pair difference timestreams. 




50 100 150 200 250 300 350 

seconds 

FIG. 3. — Timestream data for a sample scan set. The top panel shows the 
azimuth angle, and the middle panels the pair sum and difference detector 
timestreams with relative gain calibration applied; red/orange colors are the 
100 GHz pairs, while blue/green colors are the 150 GHz. The bottom panel 
shows the dark channels. In all cases an approximate scaling to temperature 
units has been applied. For the purposes of this illustration the timestreams 
have been heavily low pass filtered (to < 1 Hz). 

Figure |4] shows the mean power spectral densities (PSDs) 
of the detector timestreams. To make this plot spectra were 
computed for each channel and half-scan and then averaged 
over a day of data. We see substantial 1/f noise in the pair 
sum data. The 1 // knee shifts up in bad weather and down 
in good — the plot is for a day of intermediate quality. For 
our scanning speed (0.25°/sec in azimuth) and observing el- 
evation (~ 50°) the conversion from timestream frequency / 



to multipole on the sky is £ ~ 2000/, giving a "science band" 
of 0.1 < / < 1 Hz. Even on the worst days used in this anal- 
ysis the pair difference spectra remain close to white within 
this range. See the Instrument Paper for further details of the 
sensitivity. Note the lack of narrow line noise within the sci- 
ence band — some channels show microphonic lines at much 
higher frequencies. The roll-off above 5 Hz is imposed in the 
initial processing as part of the deconvolution procedure. Also 
note the uniform noise properties of the detectors within each 
frequency group. (The "roll up" towards higher frequencies 
observed in some pairs is due to the deconvolution of excep- 
tionally long time constant detectors.) 

3.4. Long term gain equalization via calibration source 

The QUaD telescope was equipped with a battery powered, 
remote controlled calibration source mounted behind the sec- 
ondary mirror, inside the foam cone. Before each half hour 
of observations the 45° flip mirror was commanded down and 
the polarizer grid rotated several times injecting a sinusoidally 
modulated signal into the detector timestreams — see the In- 
strument Paper for details. The low level analysis measures 
the modulation amplitude for each channel. Figure [5] shows 
the time series for one channel over the whole season and also 
the volts per airmass as measured by the elevation nods. There 
is a clear anti-correlation — as the atmosphere becomes more 
opaque the atmospheric loading goes up, suppressing the de- 
tector gains. (The source temperature is monitored and shows 
no correlation with the external temperature.) For the good 
weather data used in this analysis the suppression is < 10%. 
We make a linear regression of the mean calibration source 
timeseries against the mean elevation nod timeseries within 
each frequency band and use this to apply a correction. After 
application of this correction we believe the absolute gain of 
the QUaD system to be stable at the few percent level over the 
entire 2006 and 2007 seasons. 

3.5. Field differencing 

If we do not field difference ground pickup produces ob- 
vious artifacts in the final co-added maps. In the timestream 
this pickup is normally not visible above the detector and at- 
mospheric noise as we see in Figure [3] However on the very 
worst days, the very worst pairs can show pickup equivalent 
to as much as 5 mK CMB across a half-scan in the pair dif- 
ferenced data. We are confident that this is ground signal be- 
cause it is not present in the dark channels, is fixed in azimuth 
angle, is worst in the bottom row pixels, and correlates with 
the amount of snow present on the ground shield — see the 
Instrument Paper for further details. 

We note that the ground signal does not show structure 
smaller than a few degrees which is expected since it is near 
field pickup. It may be possible to remove the ground signal 
sufficiently well using a template based approach, and we are 
investigating this option. However for the analysis presented 
here we have applied simple lead-trail differencing through- 
out. In doing so we make the assumption that the ground sig- 
nal, which clearly does change from day to day, is stable to the 
relevant level of accuracy over the half hour timescale. In the 
field differenced maps no artifacts are visible. Ultimately the 
jackknife tests presented in Section [7] are the most sensitive 
test for residual ground contamination. 

The start/end points of each half-scan are carefully tuned 
to give a best fit to the ideal linear scan motion. We then 
apply field differencing point-by-point in the timestream by 
subtracting from each half-scan its partner occurring half an 
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FIG. 4. — Left: Timestream PSDs averaged over a single day of data. Pair sum and pair difference spectra are shown; red/orange colors are the 100 GHz pairs, 
while blue/green colors are the 150 GHz. The range of frequencies corres pond ing to 200 < £ < 2000 is enclosed within the dotted lines. See the text for further 
details. Right: The same thing for simulated noise timestream (see Section [5~2t . 
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FIG. 5. — Top: the atmospheric emission measured by elevation nods in volts per airmass for a sample detector. Middle: the sinusoidal modulation amplitude 
of the calibration source as observed by the same channel. Bottom: the calibration source amplitude after application of the loading gain suppression correction 
(see text). Histograms over all channels of the percentage fluctuation in the pre- and post-corrected timeseries are shown to the right. 
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hour later. The delta RA between the differenced points has a 
standard deviation of 0.007°. 

Note that for a Gaussian random field, such as the CMB, 
the power spectrum of a difference field (in the limit of no 
correlations) is twice that of an un-differenced field, with an 
associated increase in the sample variance due to the reduced 
effective sky area. The QUAD field differencing is explicitly 
modeled in our analysis pipeline. 

4. FROM TIMESTREAM TO MAPS 

To make polarization maps three more ingredients are re- 
quired - knowledge of the overall pointing of the telescope, 
of the relative pointing of each of the detectors, and of the 
polarization angles and efficiencies of the detectors. 

4.1. Overall telescope pointing 

The QUaD mount used a nine parameter online pointing 
model derived from optical and radio observations as de- 
scribed in the Instrument Paper. During special radio pointing 
runs this was shown to have an absolute accuracy over the 
hemisphere of ~ 0.5' rms. 

In addition a pointing check was performed on RCW38 be- 
fore and after each eight hour block of CMB observations. 
These also indicate a ~ 0.5' rms wander in the absolute point- 
ing. Attempts were made to use the measured offsets to 
make an offline pointing correction but it was not possible 
to demonstrate any clear improvement, although a few days 
were rejected due to abnormally large pointing errors. Note 
that the effect of a pointing wander of the observed magnitude 
is included in our simulations below. 

4.2. Detector offset angles 

On eleven days throughout the 2006 and 2007 seasons full 
day observations were conducted of RCW38 at three line-of- 
sight orientations ("deck" angles). For each day the data from 
each channel were fit to the six parameter model of an el- 
liptical Gaussian beam with free centroid positions, orienta- 
tion angle and widths. We did the fits both in the timestream, 
and in maps, yielding equivalent results. The scatter in the 
centroid positions for a given detector over the set of days is 
~ 0.5' consistent with the overall pointing wander of ~ 0.5' 
rms discussed above. There is hence no evidence for sys- 
tematic changes in the detector offset angles over time. We 
therefore take the detector offset angles as the mean over the 
set of observed values in the RCW38 runs and estimate their 
uncertainty as a negligible ~ 0.15'. 

4.3. Detector polarization angle and efficiency 

Our best measurements of the polarization angles of the de- 
tectors come from in situ measurements of an external source 
as described in the Instrument Paper. A chopped thermal 
source was placed behind a polarizing grid and observed with 
the telescope at many rotation angles, the signal from each 
PSB tracing out a sinusoid. The phase of these sinusoids gives 
the detector polarization angle, the fitted values agreeing with 
the design values with a scatter of around one degree rms. 
Since this is compatible with the estimated measurement un- 
certainty we have used the design values in this analysis when 
constructing maps. 

The degree to which the sinusoid mentioned above fails to 
reach zero represents the response of a measurement channel 
to anti-aligned radiation. This ratio of minimum to maximum 
response is conventionally referred to as the cross polar leak- 
age e. Our measured values of e have a mean of 0.08 with 



an rms scatter of 0.015. In this analysis we have assumed the 
mean value to apply to all channels when construc ting maps, 
and included scatter in the simulations (see Section lBTTT i. Note 
that for an experiment of this type cross polar leakage does 
not imply T to pol. leakage — it is simply a small loss of ef- 
ficiency, which must be corrected by an additional calibration 
factor applied only to the pair difference data. The effect of 
systematic error on e is discussed in Section [T0.6l below. 

4.4. Map Making 

To make maps we perform the following operations: sum 
and difference the detector timestreams for each pair, remove 
a third order polynomial across each 30 second half-scan, and 
bin into a grid of pixels weighting by the inverse variance of 
each half-scan. In this analysis the pixelization is in RA and 
Dec, and the pixels are 0.02° square. The polynomial sub- 
traction removes the bulk of the atmospheric 1 // noise allow- 
ing the simplicity of "naive" map making without incurring a 
large noise penalty. 

For the pair sum data only the signal times the weights, and 
the weights themselves, must be accumulated to form the tem- 
perature map. For the pair difference the product of the data 
and the sine and cosine of the detector angle as projected on 
the sky are accumulated. Then a 2x2 matrix inversion is per- 
formed for each pixel to produce Q and U maps. We em- 
phasize that for this matrix to be invertible any given pixel 
needs to have been measured at only two distinct grid angles 
— in the presence of noise, angles separated by 45° are opti- 
mal. The only gain from having a more uniform distribution 
of angles is in averaging down a limited number of systematic 
effects (such as beam size or pointing mismatch). In practice 
since we have detector pairs of two orientation angles, and 
observe at two deck angles, pixels in the central region of our 
maps have been measured at four angles. 

4.5. Absolute Calibration 

At this point we have T, Q, and U maps at each frequency in 
raw detector units (volts). To scale these into /iK we perform 
a correlation analysis versus two noise independent tempera- 
ture maps from the 2003 flight of the Boomerang experiment 
(B03) — which in turn have b een calibrated aga inst WMAP 
with a 2% stated uncertainty dMasi et all 12006). The B03 
maps are first passed through the QUaD simulation pipeline 
(described in Section |5TI below). The raw QUaD and "B03 
as seen by QUaD" maps are t hen b oth apodized by the QUaD 
sensitivity mask (see Section 16. U . Fourier transformed, and 
cross spectra taken between the QUaD map and one B03 map, 
and between the two B03 maps. For each bandpower b we 
then calculate our absolute calibration factor as, 

WQ.fi (m R m c ) 

ab= t r (2) 

W C fi {m R m Q ) 

where w is the Fourier transform of the beam, m are the modes 
of the Fourier transform of the apodized map, R, C and Q refer 
to the "reference", "calibration" and QUaD maps respectively, 
and the mean is taken over the modes in a given annulus of the 
Fourier plane. 

This calibration factor should ideally be a constant value 
for each bandpower (multipole range). In practice we find 
that it is for 200 < £ < 800 where the B03 beam correction 
(and therefore uncertainty on that correction) is modest, and 
we therefore take the average value across that range. We 



Second results from QUaD 



7 



have also performed the same operation using various com- 
binations of the WMAP Q, V and W band maps. Due to the 
much larger beams of WMAP, the beam corrections are very 
large, and the signal to noise low. However based on these 
results, and the point to point scatter in the B03 analysis, we 
estimate a 5% uncertainty in our primary B03 derived abso- 
lute calibration factors. 

4.6. Map Results 

Figure |6] shows the 100 and 150 GHz T, Q and U maps 
as generated using the process described in Section 14.41 The 
signal to noise in the T maps is extremely high as is evident 
from the excellent spacial correlation of the pattern between 
the two frequencies. We deliberately plot the Q and U maps 
on the same color scale to visually emphasize how small the 
polarization of the CMB is compared to the degree scale struc- 
ture in T. Three discrete sources are weakly detected in the 
temperature maps, but are not detected in polarization. 

4.7. Jackknife Maps 

The CMB polarization signal is extremely small necessi- 
tating extreme attention to detail. To probe for the presence 
of contaminating signal that does not originate on the sky we 
perform a set of jackknife tests. For each test we split the 
timestream data into two approximately equal subsets which 
should contain (nearly) the same sky signal, but which might 
contain different false signal. We then generate maps for each 
data subset and difference them to produce jackknife maps: 

M J = (M l -M 2 )/2 (3) 
y y = (y 1 + y 2 )/4 (4) 

where M J is the jackknife map, M l and M 2 are the split maps 
and V are the corresponding variance maps. Note that we 
form jackknife maps with this normalization by analogy to 
the non-jackknife map which is effectively (M l +M 2 )/2. 

Plots of the jackknife maps are not presented here as they 
are not particularly informative — they look like noise. How- 
ever Section 16.31 includes plots of the Fourier transform of 
some sample jackknife maps, and Section [7] probes for any 
hint of structure in these maps which is inconsistent with 
noise. 

5. GENERATION OF SIMULATED TIMESTREAM 

Power spectra of maps such as those shown in Figure [6] 
need to be corrected for two effects: the "noise bias" (prin- 
cipally affecting auto spectra), and the suppression of power 
due to the effects of timestream filtering and finite beam size. 
In addition we need to estimate the size of the final band- 
power uncertainties due to sample and noise vari ance. We do 
this b roadly following the MASTER technique dHivon et al.l 
2002), which requires accurate timestream level simulations 
of signal and noise. 

5.1. Signal simulations 

To construct simulations of signal timestream we start with 
LCD M power spectra generated u sing the CMBFAST pro- 
gram (Zaldarriaga&Seljak 2000) using the WMAP3 cos- 
mological paramete rs given unde r the h eading "Three Year 
Mean" in Table 2 of Sper gel et al.l (2007), and feed these into 
the "synfast" generator (part of the HEALPix package 1 ) to 

1 See http://healpix.jpl.nasa.gov/index.shtml and Gorski et alj ^20051) 



yield curved sky maps of T, Q and U at a resolution of 0.4' 
(NSIDEof 8192). 

We then read each day of data in turn, and loop over de- 
tectors. For each we calculate the sky map Md which would 
be seen by a detector of the given angle and polarization ef- 
ficiency by combining the T, Q and U input maps according 
to 

M d = i (m t + ^- (M e cos26»+M[/sin26»)^ (5) 

where e is the cross polar leakage and 9 is the detector po- 
larization angle. This ideal sky map is then convolved with 
an elliptical Gaussian smoothing kernel with the appropriate 
parameters to simulate the effect of the beam. Finally we in- 
terpolate off the smoothed map along the pointing trajectory 
for the given detector, computed as the observed telescope 
pointing direction plus the detector offset angle. 

In Section |4~31 we mentioned that fixed nominal values are 
used for the detector polarization angles and efficiencies when 
constructing maps. For each simulation realization we gener- 
ate values for each detector normally distributed about these 
nominal values with the measured rms scatter (1° and 0.015 
respectively). The simulated data is then re-mapped assuming 
the nominal values as usual. 

T o sim ulate the small pointing wander mentioned in Sec- 
tion 14.11 we generate Gaussian random numbers with zero 
mean and a of 0.5' as the RA and Dec offsets at the start 
and end of each eight hour block. We then linearly interpo- 
late these offsets to each time-step and add them in to the 
observed pointing trajectory. Although in reality the wander 
has a more complex behavior given the observation strategy, 
and averaged over hundreds of days of data, this will lead to 
a broadening of the effective beam width of the correct (very 
small) amount. 

The detector offset angles are taken as fixed at our best esti- 
mate values — we ha ve no evidence for time variation as men- 
tioned in Section l4~2l above. Note that the measured beam cen- 
troid positions show repeatable offsets between the two halves 
of each detector pair with rms magnitude of ~ 0. 1 '. When we 
construct maps using pair sum and difference data we use the 
mean for each pair. However when we sample from the sky 
to generate simulated timestream we use the measured indi- 
vidual detector values. Hence any T to pol. mixing which 
occurs due to beam centroid mismatch is fully included in the 
simulations — and found to be negligible — see Section [T0.4l 

Note that generation of simulated signal timestreams takes 
place before the field differencing operation — effects such 
as actual signal correlations and mismatch of pointing coordi- 
nates between lead and trail fields will therefore be included. 

5.1.1. Measurement of the beam widths 

To convolve the ideal sky map with the beam we need to 
know the beam shape for each detector. We have obtained this 
from nine days of special observations of the bright quasar 
PKS0537-441, which is a point source at th e angular resolu- 
tion of our experiment dFev & Charlotl2000l) . The mean of the 
measured major and minor FWHM's is 5.0' at 100 GHz and 
3.5' at 150 GHz with uncertainties of w 2%. There is evidence 
for a small degree of variation in width between the detectors 
of a given frequency band, and for < 10% elongation (mis- 
match between the major and minor widths), but for any given 
detector these are sub-dominant to the measurement uncer- 
tainty. As when determining the detector offset angles from 
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FIG. 6. — QUaD 100 and 150 GHz T, Q and U field difference maps. The color scale is ±200 fiK in all cases, and the maps have been smoothed with a 5' 
FWHM Gaussian kernel. (Black circles indicate discrete sources which are removed in the power spectrum analysis.) 
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measurements of RCW38 (see Section 14.21 above) we have 
performed the elliptical Gaussian fits in both the timestream 
and in maps yielding nearly identical beam widths and angles. 
For further details see the I nstrument Paper and our Optics Pa- 
per (lO'SullivanetalfeOOSl) . 

The effective beam width in our CMB field co-added maps 
is measured with low signal to noise by the quasar PKS0524- 
485 (the brightest of the circled sources in Figure [6). This 
source appears as a Gaussian peak with width consistent with 
the single day observations of PKS0537-441. 

5.2. Noise simulations 

As shown in Figures [3] and [4] the detector timestream is 
dominated by heavily correlated low frequency atmospheric 
noise. The goal of the timestream noise generator is to repro- 
duce simulated half-scans which are indistinguishable from 
the real under a battery of tests in both the time and frequency 
domains. 

To achieve this we find it necessary to measure and re- 
generate longer pieces of timestream, and then cut them down 
to the half-scans which are actually used when constructing 
the maps. First the complete time period spanning each five 
scan set is Fourier transformed (including turn arounds). Then 
for each of a set of logarithmically spaced frequency bins 
we take the covariance matrix of the Fourier modes between 
all channels. This matrix is Cholesky decomposed and used 
to mix uncorrelated random numbers to re-generate the ob- 
served degree of covariance. This process is repeated for each 
frequency bin and then the resulting sets of Fourier modes 
are inverse transformed to yield simulated timestream. Pos- 
sible correlations between the real and imaginary parts of 
each Fourier mode are preserved by using complex covariance 
and Cholesky matrices, but this process assumes the Fourier 
modes are correlated only between channels, and not between 
modes. Breaking the simulated timestream down into half- 
scans reintroduces such correlations and we find it has equiv- 
alent temporal and spectral characteristics to the real data un- 
der a variety of statistical tests. 

As an example the right part of Figure H] shows the power 
spectra of the resulting simulated timestream — comparing 
to the real PSDs in the left part of the figure we see an ex- 
cellent match. Note that the simulation produces non pair 
differenced timestream which has then been differenced to 
generate the spectra shown in Figure [4] — it is clear that the 
channel-channel correlations are being reproduced to high ac- 
curacy. The narrow lines in the real spectra are re-generated as 
broader boxcar features due to the frequency bin width used 
in the simulation process — this is irrelevant since these fre- 
quencies are far above the science band. 

The signal to noise ratio in the timestream data is suffi- 
ciently low that we do not need to subtract the signal con- 
tribution before taking the noise spectra. (The PSD of signal 
only LCDM timestream is more than 2 order of magnitude 
below the noise at all frequencies for the pair sum, and 4 or- 
ders of magnitude for the pair difference.) To confirm that the 
noise modeling operates correctly we have performed a "sim 
the sim" study where the output of a single signal plus noise 
timestream simulation is used as the input to a complete simu- 
lation set, passing all the jackknife tests described in Section|7] 
as expected. 

5.3. Simulated maps 

We generate many realizations of signal and noise 
timestream and produce the full set of un-split and jackknife 



split maps from these using the same code used to make maps 
of the real timestream in Section 14.41 above. For each real- 
ization we then add the signal and noise maps to form signal 
plus noise maps — since the map making process is linear this 
is equivalent to adding the timestreams, and computationally 
more efficient. The generation of simulation realizations is 
computationally costly and hence their number is relatively 
small (500 in this analysis). 

6. FROM MAPS TO POWER SPECTRA 

This section will describe the various steps which take us 
from the real and simulated maps to angular power spectra. 

6.1. The apodization mask 

In addition to the signal maps Mj q u the process described 
in Section l4~4l also produces maps of estimated variance Vt,q,v 
based on the variance of each co-added half-scan, assuming 
that the noise is white. Empirically these are found to be close 
to correct for the polarization maps and a 10% (20%) underes- 
timate for the 100 GHz (150 GHz) T maps, as expected since 
the pair sum data contains significant 1/f noise. We take the 
inverse of these variance maps as the apodization masks: 

Ax = ^ (6) 
Vx 

However due to the partially overlapping "tiles" of cover- 
age resulting from our observation strategy and the differ- 
ing detector offset angles there are sharp steps in the masks 
(of ~ 10% magnitude). If these steps are not smoothed out 
then the product of the map and mask (the quantity which is 
about to get Fourier transformed) will also have sharp steps. 
Such steps correspond to a mixing of power from large to 
small scales which is particularly undesirable in the TT spec- 
tra where even a small fractional contribution from the very 
large low multipole Q values would dominate over the intrin- 
sic power in the damping tail region of the spectrum 2 . To mit- 
igate this effect we convolve the inverse variance mask with a 
Gaussian shaped smoothing kernel of FWHM 0.5°: 

A' X =A X *G (7) 

We emphasize that an arbitrary apodization mask can be 
used without biasing the results and, in fact, the inverse vari- 
ance mask is only optimal in the limit of low signal to noise 
which is not the case for the T maps. In principle an optimal 
mask could be generated for each spectrum, and in fact each 
bandpower within each spectrum, based on its signal to noise, 
but we have not pursued this complication in this analysis. 

After smoothing we inject Gaussian shaped "divots" with 
FWHM of 0.5° into the apodization masks at the locations of 
three discrete sources circled in Figure[6]to null out any effect 
they might otherwise have on the results: 

where x and y are the pixel coordinates, x„ , and y 0> i are the 
source locations, and the loop i is over the set of sources to be 
nulled. Note that although the sources are visible only in the 
T maps we apply the same masking in polarization also. 

2 Note that the relevant ratio is in Q, not £(£+ 1)Q which is the quantity 
conventionally plotted — for the TT spectrum under LCDM C200/C2000 ~ 
2500 
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The final apodization masks imply an effective sky area for 
this analysis of 25 square degrees. 

6.2. Fourier Transform and Power Spectra 

We next make the fiat sky approximations and take the two 
dimensional discrete Fourier transform of the product of the 
map and apodization mask: 

m x =cVY(M x A'i) (9) 

where c is a normalization constant which acts to make the 
apodization and Fourier transform operations total power pre- 
serving. 

In the Fourier domain the transform from Q,U to E,B is 
simply 

mE = +mQCOs2cf)+musiii2(f> (10) 
m B = -m e sin2(/)+OTj/Cos20 (11) 

where <f> is the polar angle of each mode m with respect to the 
Fourier plane origin. 

We then take products of the (complex) modes within and 
between each set 

pxY = m x mY (12) 

where X and Y can be T, Q or U at either 100 or 150 GHz. 

Each set of mode products pxy is then multiplied by d = 
£(l+l)/2w where the multipole £ = 2iru, and u is the Fourier 
conjugate variable to angular distance from the map center. 
Finally we take the binned angular power spectra 

bxr,i= (dpxY,i) (13) 

where the mean is taken over the modes in each annulus i of 
the Fourier plane. The multiplication of the map by the mask 
in real space corresponds to a convolution in Fourier space 
— hence the Fourier modes mx are correlated and so are the 
bandpower estimates bxY- For this analysis we set the bin 
spacing to A£ = 81 which results in correlation of ~ 20%. No 
importance should be attached to the exact value 81 — it was 
an arbitrary choice made early on and remained unchanged 
throughout the development of the analysis. 

The above description is a slight simplification of the ac- 
tual procedure used. In reality for each cross spectrum XY 
we enforce a common apodization mask A'^ Y = ^jA'^Ay be- 
fore taking the Fourier transforms and proceeding to power 
spectra. 

6.3. Fourier Plane Masking 

The strong theoretical expectation is that the CMB is 
isotropic on the sky and therefore also in the two dimen- 
sional Fourier plane. The noise however is certainly not — 
atmospheric variation injects strong low frequency fluctua- 
tions (1// noise) along the scan direction and hence in a band 
around the perpendicular direction in the Fourier domain. We 
are perfectly at liberty to down weight or excise portions of 
the Fourier plane when we take power spectra — no bias will 
result so long as the weight/cut criteria are independent of the 
mode values of the real maps. 

Figure|7]shows the 150 GHz T and Q Fourier modes for the 
deck jackknife map. The third order polynomial filter which 
has been applied to each half-scan removes power along the 
scan direction producing the dark vertical bands down the 



middle of the plots. Atmospheric 1/f which survives the fil- 
tering appears as a brighter band either side of the dark band 
in the T plot, but due to the unpolarized nature of the atmo- 
sphere, does not appear in the Q plot (see Section l5T2l above). 
Because the atmospheric noise is highly correlated amongst 
adjacent detector pairs in the array we in fact see a complex 
structured pattern of non-uniform noise in the T plot — our 
noise simulations fully reproduce this pattern. 




-2000 -1000 1000 2000 

27CU 



FIG. 7. — Fourier transforms of the apodized 150 GHz T and Q deck jack- 
knife maps. The real part of the square of the Fourier modes is shown with 
a linear color stretch. The red circles indicate the annuli within which means 
are taken to generate bandpower values. Note that the horizontal direction in 
this plot corresponds to the scan direction — see text for further details. 

Section loTTl mentioned the need to smooth the inverse vari- 
ance derived apodization mask to reduce up-mixing of power 
from small to large multipoles. There is an additional simi- 
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lar, but smaller, effect. Due to the detector offset angles each 
has observed a slightly different rectangular patch of the CMB 
sky. We polynomial filter each detector's timestream for each 
half-scan to reduce atmospheric noise, and so the CMB modes 
removed from each "tile" are slightly different. Hence when 
we co-add the data from all detectors small step artifacts are 
present in the overall map at the tile overlap boundaries. Ex- 
amining the non-jackknife Fourier T modes using plots simi- 
lar to those shown in Figure|7]we see up-mixed power becom- 
ing dominant over the intrinsic for multipoles at I > 1800 in 
a narrow band around the scan direction axis — as expected, 
since the artifacts are vertical step edges in the maps. This 
spurious power is also seen in the signal only simulations con- 
firming its origin is up-mixing. Note that this effect is only 
significant in T due to the much larger ratio between low and 
high multipoles than in the polarization spectra. However we 
excise these modes when generating all power spectra result- 
ing in a trivial increase in the uncertainty of our highest band- 
powers (< 3%). This effect will be more important in a future 
analysis extending to £ = 3000. 

One could weight the Fourier modes when combining to 
form power spectra based on their signal to noise ratio as mea- 
sured in the signal and noise simulations. This would be (at 
least a partial) substitute for the polynomial half-scan filtering 

— one either removes the atmospheric 1 // in the timestream 
or in the two dimensional Fourier plane. Looking at Figure|7] 
such Fourier plane weighting would clearly be beneficial for 
the TT spectra. Although we have experimented with such 
weights and cuts we have not implemented them in the analy- 
sis presented here — the benefit for the essentially white noise 
polarization spectra is close to zero. 

6.4. Raw spectra compared to simulations 

By taking the mean of the masked Fourier mode products 
within each annulus we generate raw power spectra of the real 
and simulated maps. Figure [8] compares the 150 GHz EE 
spectra and is in a sense the fundamental result of our analysis 

— is the observed spectrum consistent with being a realization 
of LCDM plus noise? To test that hypothesis we could sim- 
ply perform \ 2 tests at this point. However, for presentation 
purposes, and to allow our spectra to be used for cosmologi- 
cal parameter analyses, we proceed to noise and filter correct 
them as follows. 

6.5. Removing the noise bias 

At this point in the analysis we have sets of bandpower val- 
ues calculated using Equation [13] for each of the auto and 
cross spectra. We have these for the real spectra, and each 
realization of the signal only, noise only, and signal plus noise 
simulations. We will denote these Txy, Sxy, Rxy and sn^y 
respectively. 

To remove the noise bias we take the mean of the ensemble 
of noise only simulations and subtract from the real spectra 

r' XY = r XY -nxY (14) 

In Figure|8]the blue curve is subtracted from the black points. 
Correlated noise, presumably from the atm osph ere, is fully 
simulated by the process described in Section l5T2l above. Such 
noise can result in the mean noise cross spectra — where cross 
spectrum means either within, or across frequency bands — 
being non zero. In practice the levels are very small compared 
to the auto spectra but for completeness we subtract them any- 
way. 
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FIG. 8. — Comparison of real and simulated raw 150 GHz EE spectra. 
The cyan curves are the ensemble of signal plus noise simulations, the black 
points are the observed bandpower values, and the blue curve is the mean of 
the ensemble of noise only simulations. 

6.6. Band power window functions 

As mentioned in Section |oT2l multiplication by the apodiza- 
tion mask in image space corresponds to convolution in the 
Fourier plane by the Fourier transform of that mask. In Fig- 
ure [7] we see that adjacent Fourier modes (pixels) are highly 
correlated, and that the resulting annular mean bandpower es- 
timates will therefore also be correlated. We thus need to 
know how much each multipole on the sky contributes to each 
experimental bandpower — the so cal led "band power win- 
dow functions" (BPWF's) dKnox| [l999). This could be deter- 
mined by running sets of signal only simulations, each with 
non-zero input power only in a narrow band of multipoles. 
However the computational cost of doing this in practice is 
too high, and we hence use an alternate, much faster method. 

We take each narrow annulus in the Fourier plane in turn 
and convolve this with the Fourier transform of the mask. For 
Q and U we generate annuli corresponding to pure E or B 
modes. Taking the power spectrum of the convolved annuli 
measures the response of each experimental bandpower to sky 
power in the given annulus. Placing these spectra into the 
rows of a matrix, the columns are then the desired BPWF's 
(interpolating to every multipole). For BB we have two func- 
tions per bandpower representing the response to true BB sky 
power, and also the response to EE. In practice the mixing 
is ~ 10% in our lowest bandpower and falls rapidly with in- 
creasing I (see Section flO.lt . Figure [9] shows some example 
BPWF's. 

6.7. Determining the filter/beam suppression factor 

Using the BPWF's described above we can calculate the 
expectation values e^y of the observed bandpowers given the 
input LCDM power spectra used to generate the signal simu- 
lations. We then take the ratio of the mean of the signal only 
simulated spectra to these expectation values to empirically 
determine the suppression of power which has occurred due 
to the convolution of the sky by the telescope beam pattern, 
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FIG. 9. — The bandpower window functions for the 150 GHz BB spec- 
trum. Each blue curve shows the relative response to BB sky power, while 
the corresponding red curves show the response to EE. 

and the polynomial filtering of the timestream 



— this is the factor by which the real spectra must be divided 
to yield an un-biased estimate of the true sky power. Note that 
this process will automatically include any simulated suppres- 
sion effect (including the so-called "pixel window functions"). 
Figure [10] illustrates this step — the curve in the lower panel 
approaches the "beam window function" Wi at higher £. 
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FIG. 10. — Upper: Input EE spectrum (red line), calculated expectation 
values (red crosses), and mean simulated signal only spectrum (magenta cir- 
cles). Lower: The ratio of the two sets of points in the upper panel — the 
filter/beam suppression factor. (This plot is for the 150 GHz spectrum.) 



We now divide the spectra by their respective suppression 
factors 

r' 

t xy = j— (16) 

In practice we apply the fjr correction to the TT spectra and 
the (ee correction to the EE, BB and EB spectra. For the TE, 
TB spectra we use the geometric mean of fjj and fg£. 

6.8. Power spectrum results 

We also apply the noise de-bias and filter/beam correction 
operations of Equations [14] and [16] to each signal plus noise 
realization. Their fluctuation then provides an estimate of the 
uncertainty of the real bandpower values. This uncertainty 
will only be correct in as much as the theory spectrum used 
as input to the signal simulations matches reality — if there is 
any significant disagreement one should iterate the enti re pro- 
cess until it converges. As we will see in Section |9^2l our re- 
sults are in fact perfectly compatible with the WMAP3 based 
model used as input to the signal simulations (see Section[5j} 
so there is no need to iterate. 

Figure [Til shows the full set of 21 signal spectra: within 
each frequency band we have TT, TE, EE, BB, TB and EB, 
while across frequency bands we can form TiooTiso, 7ioo£i50> 
£ioo7i50> E\oaE\sQ, B\qqB\so, ^ioo^iso, ^100^150, ^ioo^iso and 
2?iooi?i50- We see that the 150 GHz spectra have better sen- 
sitivity than the 100 GHz — this is partly due to the larger 
number of detectors at the higher frequency, but mostly due 
to the smaller beam size (and hence smaller beam correction). 

6.9. Alternate analysis 

The main analysis in this paper come s from an evo lution 
of "pipeline 2" in our previous paper lAde et all d2008l) . A 
second pipeline exists, the low level parts of which are de- 
rived from our previous "pipeline 1", but which, for the re- 
sults presented here, is now also performing flat sky power 
spectrum estimation. This pipeline currently has noise mod- 
eling that is so mewhat less sophisticated than that described 
in Section 15.21 The two pipelines were independently writ- 
ten and share no code. Although the algorithms implemented 
are intended to be basically the same they are sure to differ in 
some (hopefully unimportant) details. The fact that the final 
results agree very closely therefore adds considerable addi- 
tional confidence. Figure [12] compares the 150 GHz spectra 
from the two pipeline. 

7. JACKKNIFE TESTS 

To probe the results shown in Figure QT] for systematic 
contamination we conduct a battery of jackknife tests. As 
mentioned above split dataset maps of the real data and each 
simulation realization are generated. Each split divides the 
timestream data into two approximately equal subsets which 
should contain (nearly) identical sky signal, but which might 
contain different contaminating signal. In this analysis we 
consider four such data splits which we call the deck angle 
jackknife, scan direction jackknife, split season jackknife and 
focal plane jackknife. Each of these is described in more de- 
tail below. 

For each jackknife the split maps are subtracted from one 
another, divided by two, and we then proceed to power spec- 
trum estimation as usual. For the real spectra we then calcu- 
late the x 2 versus the null model 
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FIG. 1 1 . — The full set of QUaD power spectrum results. The blue points are the signal spectra, the magenta the deck jackknife (see text for definition), and 
the red curves a conventional LCDM model. The horizontal and vertical axes are multipole £ and £(£+ 1)C; /2tt ([iK 2 ) respectively. Each bandpower is ~ 20% 
correlated with its neighbors and, where signal dominated, there are strong correlations between bandpowers across frequency bands. The error bars are the 
standard deviation of the simulated signal plus noise bandpower values, and the two sets of points have been offset by ±15 in £ from their nominal values for 
clarity. Note that absolute calibration and beam size uncertainty are not included in the error bars. The jackknife spectra are consistent with null — see TablefT] 
and the text. 
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500 1000 1500 2000 500 1000 1500 2000 

FIG. 12. — Comparison of 150 GHz spectra for the main (blue points) 
and alternate (magenta points) pipelines. The horizontal and vertical axes are 
multipole t and l(£ + 1)C/ /2tt (/iK 2 ) respectively, and the two sets of points 
have been offset by ±15 in i from their nominal values for clarity. 



x 2 = r"C' l r'" (17) 

where C is the bandpower covariance matrix as estimated 
from the ensemble of signal plus noise simulations. We also 
calculate \ 2 for each signal plus noise simulation. The gener- 
ation of simulation realizations is computationally costly and 
hence their number is relatively small (500 in this analysis). 
This has two implications. Firstly, since bandpower correla- 
tions beyond nearest neighbor are sufficiently weak as to be 
lost in the measurement noise, we set all but the main and 
first two off diagonals of the bandpower covariance matrix to 
zero. Secondly when calculating the \ 2 values for the signal 
plus noise realizations extreme fluctuations have the opportu- 
nity to partially self-compensate by injecting extra covariance 
into the matrix which they will be measured against, biasing 
the resulting \ 2 values low. To get around this we re-calculate 
the covariance matrix excluding each realization in turn be- 
fore calculating the \ 2 f° r that realization. 

We find that the simulated jackknife \ 2 distributions for 
TT do not follow the analytical expectation. Examining Fig- 
ure [13] which shows the bandpower deviations of the deck 
jack spectra, we see that there are two reasons why. Firstly 
the simulated bandpower distributions are significantly non- 
Gaussian — the colored lines do not sit at -2 through +2. 
Secondly, and more importantly, there is significant predicted 
imperfect cancellation — the simulated distribution deviates 
strongly from a median of zero at large angular scales. The 
reason for this latter effect was mentioned in Section [531 — 
in any split where the sky coverage "tiling" for the two sub- 
sets is non-identical the interaction of the true sky brightness 
distribution and the half-scan polynomial filtering generates 



TABLE 1 

Jackknife PTE values from \ 2 tests 



Jackknife 


100 GHz 


150 GHz 


Cross 


Alt. Cross 


Deck angle 










TT 


0. 102 


0.072 


0.030 




TE 


0.066 


0.040 


0. 174 


0.586 


EE 


0.874 


0.624 


0.210 




BB 


0.990 


0.316 


0.946 




1 D 


U. jUO 


f\ A 1 A 


U. jzU 


U.oiz 


EB 


0.736 


0.210 


0.952 


0.180 


Scan direction 










TT 


0.650 


0.238 


0.982 




TE 


0.316 


0.444 


0.870 


0.718 


EE 


0.656 


0.490 


0.372 




BB 


0.882 


0.156 


0.982 




TB 


0.274 


0.430 


0.154 


0.254 


EB 


0.540 


0.332 


0.826 


0.418 


Split season 










TT 


0.202 


0.020 


0.408 




TE 


0.962 


0.296 


0.300 


0.472 


EE 


0.116 


0.038 


0.066 




RR 


U.oVD 


U.UjO 


n i on 

U. 1UU 




TB 


0.946 


0.494 


0.896 


0.166 


EB 


0.860 


0.612 


0.972 


0.890 


Focal plane 










TT 


0.352 


0.472 


0.232 




TE 


0.858 


0.872 


0.448 


0.936 


EE 


0.654 


0.174 


0.944 




BB 


0.696 


0.642 


0.708 




TB 


0.888 


0.422 


0.008 


0.456 


EB 


0.932 


0.640 


0.354 


0.344 



slightly different output maps. Consistent with our overall 
simulation based approach we therefore take the probability 
to exceed (PTE) the real spectra \ 2 values versus the simu- 
lated distributions, rather than the analytical. 

Table [JJ shows the full set of PTE values for all the jack- 
knifes and spectra — there are no strong indications of prob- 
lems. These values are expected to be uniformly distributed 
between zero and one, and in Figure [14] they are seen to be 
so. See the following sub-sections for detailed discussion of 
these results. 

7.1. Deck jackknife 

The deck angle jackknife is perhaps the most powerful. 
Due to the locking of our daily observations to LST the two 
halves of this split contain data taken over completely dif- 
ferent ranges of telescope azimuth. In addition, and equally 
importantly, the entire telescope is rotated by 60° around the 
line of sight between the two observation sessions. The two 
data subsets are therefore separated from one another in time, 
azimuth angle, and detector polarization angle as projected on 
the sky. It is very hard to conceive of a source of contamina- 
tion which would be common in the polarized maps for the 
two halves of this split. We see no significant indication of 
problems with the deck angle jackknife in Figures [TT] and [13] 
or in Table [JJ One might perhaps worry about the rather low 
values for 150 GHz and cross-frequency TT but we caution 
against over-interpreting such numbers — the table contains 
84 numbers so on average four numbers below 0.05 are to be 
expected. In addition looking at Figure [13] we see that the 
bandpower making the strongest contribution to the cross fre- 
quency x 2 has a very low deviation in the 150 GHz spectra. 
Note that Figure[l4]shows no obvious excess of low PTE val- 
ues. 
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FIG. 13. — Bandpower deviations of the deck jackknife spectra. The horizontal axis is multipole number. The spectra shown as magenta points in Figure [Til 
have been divided by their error bars to investigate the contribution of each bandpower to \~ ■ The green, red and blue lines show the 2.3%, 15.9%, 50%, 84.1% 
and 97.7% points of the integral distribution of the simulated signal plus noise bandpower values. 
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FIG. 14.— Distribution of the \ 2 PTE values from Table[TJ 

7.2. Scan direction jackknife 

The scan direction jackknife splits the data into the out- 
going and returning half of the scans. In terms of external 
contamination this is perhaps the easiest test to pass — only 
something very rapidly varying would cause it to fail. How- 
ever it is a stringent test for internal instrumental effects. Any 
scan synchronous false signal, caused perhaps by motion of 
the liquid cryogens in the tanks provoked by the telescope 
motion, would likely fail to cancel in this jackknife. Also any 
failure to adequately deconvolve the temporal response of the 
detector channels would cause this test to fail. Looking at 
Table Q] we see no problems. 

7.3. Split season jackknife 

For the split season jackknife we have divided the time or- 
dered list of days in half (not in fact into 2006 and 2007 — 
the split occurs in late August 2006). We then make maps 
with each set of days and difference them. This test would 
fail if there were a significant shift in the absolute calibration 
of the telescope system beyond the atmosp heric loading ef- 
fect whose correction is described in Section [3~4l Looking at 
Table Q] we see a few numbers below 0.05 but nothing highly 
significant. Examining the bandpower deviations (a plot not 
shown in this paper analogous to Figure [T"3l) we find that for 
each spectrum with a low PTE it is caused by a different band- 
power(s), giving no further hint of any problem. 

7.4. Focal plane jackknife 

This jackknife splits the detectors into the two orienta- 
tion groups which are separated by 45° — referred to as 
instrument-^ and instrument-?/. Due to co-adding across the 
two deck angles each map pixel has still been observed at 
two polarization angles allowing the construction of Q and 
U maps as usual. This test is perhaps the weakest, but might 
reveal problems with instrumental false signal in a subset of 
the pairs (although it is hard to see how that would not also 
show up in the deck angle jackknife). Looking at Table Q] we 
see no problems. Although this jackknife includes the lowest 
number in the table (0.008) one such number is not improb- 
able and we note again that the PTE distribution shown in 
Figure[l4]is consistent with uniform. 

8. FOREGROUND STUDIES 



8.1. Frequency difference maps and spectra 

We can also take the difference between our 100 GHz and 
150 GHz maps. It is important to be clear that this is not a 
jackknife in the sense of Section|7] — the true sky brightness 
distribution may in fact differ at these two frequencies. There- 
fore any failure to cancel might be due to the presence of as- 
trophysical foregrounds, as well as instrumental systematics 
or contamination. 

As described in Section l4~5l we find the absolute calibration 
scalings for our 100 and 150 GHz maps by cross correlating 
them against the same B03 maps (which are at 150 GHz). 
However, if the sky pattern differs at these two frequencies 
then this will still show up in the frequency difference maps. 
Figure Q3] shows the difference between the real 100 and 
150 GHz T maps, and the same thing for a signal plus noise 
simulation realization. 

As mentioned in Section|7l when subtracting maps with dif- 
ferent sky coverage "tiling" a small degree of mismatch is ex- 
pected due to the polynomial filtering. This is the cause of the 
vertical "step edges" observed in the upper panel of Figure[T"5l 
The simulation realizations show similar effects, although in 
the real map they do appear to be unusually strong. Note that 
we also expect to see cancellation failure at smaller angular 
scales due to the differing beam sizes for the two frequency 
bands. 

We next take power spectra of the frequency difference 
maps and compare them to simulation as was done for the 
jackknifes in Section [7] Figure [T6l shows the resulting band- 
power deviations — PTE values analogous to those in Ta- 
bleware shown in the Figure. The TT spectrum shows larger 
deviations at low multipoles than expected from the simula- 
tions, and the probability that these are caused by the differ- 
ing "tiling" effect alone is low. However the absolute value of 
these bandpowers is ss 15 [iK 2 to be compared to the 1000's 
of /iK 2 in the un-differenced map — i.e. the fractional can- 
cellation failure is very small, and completely irrelevant com- 
pared to the sample variance in the TT spectra. It is not clear 
whether the excess cancellation failure is due to instrumental 
effects or real foreground signal. 

8.2. Predicted diffuse foreground levels 

Our field was chosen to partially overlap the B03 deep field 
to allow absolute calibration against that map. They in turn 
chose the location based in part on the position of the Sun dur- 
ing their balloon flight. Although low, the foreground emis- 
sion in this region is not the lowest available on the sky. 

To obtain estimates for the expected level of foreground 
dust e mission we use the FD S model 8 extrapolation of IRAS 
maps (Finkbein er et al.l 19991) — this is shown in Figure[2] For 
synch rotron emission we use an extrapolat ion of 408 MHz 
maps ( Haslam et al.ll 198 it lFinkbeinerll20 01). We pass these 
maps through the QUaD simulation pipeline, including the 
field differencing and filtering operations. For dust the max- 
imum of the resulting TT spectra is ~ 4 /iK 2 (150 GHz, t = 
200) while the synchrotron maximum is a negligible 0.03 /iK 2 
(100 GHz, £ = 200). Although these models are possibly not 
the current best available data it is clear that such extrapola- 
tions will not give detectable levels in either temperature or 
polarization, in the presence of CMB, and given the sensitiv- 
ity of QUaD. 

8.3. Point Sources 

Three discrete sources are visible in the 100 and 150 GHz 
T maps shown in Figure|6] As described in Section loTTl these 
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FIG. 15. — Upper: The difference between the 100 and 150 GHz T maps 
shown in Figure |6]on a ±20 pK color stretch. Lower: The same thing for a 
signal plus noise simulation realization. 

are masked before taking the power spectra. Turning off the 
masking of these sources we find that the 100 GHz TT band- 
powers increase by ps 10% at the highest multipole consid- 
ered in this paper (2000) with the increase falling off to lower 
multipoles. The 150 GHz TT spectrum shows an increase 
of ps 3% at t of 2000 and a similar falloff. The EE and BB 
bandpowers show changes of <C 1%. 

The source flux distribution dN/ds is typically a power law 
with the majority of the anisotropic power being contributed 
by the brightest few sources. We therefore estimate the resid- 
ual point source contribution to be < 3% in the highest bin of 
100 GHz TT, < 1% in the highest bin of 150 GHz TT and 
negligible in all other spectra. 

8.4. Template cross correlation 

To test for the possibility of emission correlated with ther- 
mal dust, but stronger than expected on the basis of extrapola- 
tion, we have carried out a template cross correlation study. 
After passing the FDS dust maps through our pipeline the 
resulting maps were cross correlated with the corresponding 
QUaD CMB maps. Though we might expect some non-zero 
corre lation simply by a c hance alignment of large scale struc- 
ture dChiang et al.l 12007). we find that compared to simula- 



5000 
4000 
3000 
2000 
1000 


100 
50 

-50 
-100 
-150 
40 

20 



-20 
20 
10 


-10 
-20 






•••••••••••••••••••••a* 



II 





t PTE = 0.000 



■/!' PTE = 0.558 




X PTE = 0.772 




X PTE = 0.336 



X PTE = 0.374 



X PTE = 0.644 



1000 1500 



FIG. 16. — Left: Frequency difference spectra and Right: the associated 
bandpower deviations. \ 2 i s taken using the bandpower covariance matrix 
and converted to the PTE value shown on each pa nel. (For further explanation 
see the captions to the analogous FiguresfTT]and[T3]) 

tions there is no evidence of contamination. Indeed analyzing 
the cross power spectra at both 100 and 150 GHz between 
dust and CMB on a per bandpower basis, against simulations, 
reveals no problems, nor any suggestion that dust foregrounds 
are responsible for the cancellation failure seen in the lower 
bandpowers of the frequency difference TT spectrum. 

9. COMBINED SPECTRA 

We now wish to form a single combined set of spectra from 
the 100 GHz, 150 GHz and frequency-cross results presented 
above. For each bandpower of each spectrum we take the 
3x3 covariance matrix over the ensemble of signal plus noise 
simulations (4x4 for TE, TB and EB). The combination 
weights are the column (or row) sums of the inverse of this 
matrix. The improvement over the 150 GHz bandpower un- 
certainties is between zero and 30% depending on whether 
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the bandpower is signal or noise dominated. The BPWF's 
are also combined. Figure [17] shows the combined spectra 
as compared to their expectation values under LCDM. The 
plotted bandpowers, together with their covariance matrices 
and bandpower window functions, are available in numeri- 
cal form at http://quad.uchicago.edu/quad In contrast to the 
jackknife spectra, for the signal spectra we find the simulated 
bandpower distributions to be Gaussian (Cf. Section|7]i. 
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FIG. 17. — Combined QUaD power spectra shown as black points with 
error bars. The red crosses are the expectation values for each bandpower 
given the LCDM model plotted as a red line. Note the differing y-axis scales. 



9.1. Absolute Calibration and Beam Systematics 

In addition to the sample and noise variance error bars 
shown in Figure [17] there are two major sources of system- 
atic uncertainty. The first of these is the uncertainty on our 
absolute calibration against the B03 maps. As mentioned in 
Section l431 we estimate this uncertainty as 5% in temperature 
units (10% in power). To estimate this uncertainty one simply 
multiplies the bandpower expectation values by 0. 1 , and takes 
the outer product as an addition to the bandpower covariance 
matrix: 



C fl : 



(18) 



where a a = 0.1. 

The second major systematic effect is uncertainty on the 
width of the telescope beam. It is highly unlikely that the 
widths used in the simulations are significantly broader than 
the true values — they are very close to the results of physi- 
cal optics calculations (O' Sull ivan et alj|2008l) . However, it is 
conceivable that we have somehow under-estimated the point- 
ing wander and associated effects discussed in Section 15.11 
above, causing the effective overall beam width in the simu- 
lations to be narrower than that in the real maps, and the sup- 
pression factor curve plotted in Figure [10] to be higher than it 
should be. This would cause the corrected bandpower values 
to be biased low with increasing £. 

In addition we have not yet carried out exhaustive inves- 
tigations of the beam shape and measurement uncertainties. 
Therefore we very conservatively assign a beam width uncer- 
tainty of 10% for this analysis with the expectation that this 
will be improved upon in a future analysis dedicated to high 
I TT. The effective beam FWHM for the combined spectra 
is 4.1' and the fractional shift in the bandpower values which 
would result from increasing this is well approximated by 



S e = 



We 



.] = g ^+2W+l)_l 



(19) 



where 07, = 6>fwhm/ \/81n2, and S = 0.1 is the fractional beam 
error. To estimate this uncertainty one multiplies the band- 
power expectation values by the Si factors calculated at the 
band center I values, and takes the outer product as an addi- 
tion to the bandpower covariance matrix. 



(20) 



Figure [18] shows the magnitude of these uncertainties, as 
well as the sample and noise variance, for each of our com- 
bined spectra. All spectra except BB are sample variance 
dominated at lower £, and all but TT are noise dominated at 
high I. Absolute calibration uncertainty is sub-dominant to 
random uncertainty for all spectra at all I. Beam uncertainty 
becomes the dominant effect for TT at high I. 

9.2. Comparison to LCDM 

As mentioned in Section [9] for the signal spectra we find 
the simulated bandpower distributions to be Gaussian. In ad- 
dition we find the simulated \ 2 distributions to be consistent 
with the analytical expectation. Therefore, in contrast to the 
jackknife spectra, for the signal spectra we quote the analyt- 
ical probability to exceed. Figure [19] shows the bandpower 
deviations comparing the combined spectra to LCDM and the 
null model. Our results are clearly perfectly compatible with 
LCDM and crushingly incompatible with the no polarization 
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FIG. 18. — Contributions of the various sources of uncertainty for the com- 
bined bandpowers. 

hypothesis. Note that ~ 10 EE bandpowers have > 4a signif- 
icance. 

9.3. Investigation of EE peaks 

Looking at Figure [T7] we appear to see three or four of the 
expected acoustic peaks in the EE spectrum. It is interest- 
ing to ask at what significance these have been detected. One 
way to do this is to determine the \ 2 of the observed band- 
powers against a version of the LCDM model which has been 
smoothed sufficiently to remove the peaks. Figure l20l shows 
the result — the LCDM model has been convolved with a 
Gaussian with a> = 150. We see that the probability that such 
a model is correct is very low, the \ 2 PTE being 0.001. 

In LCDM the TT and EE peaks are approximately half a 
cycle out of phase with one another. If l s is the peak spacing, 
and n is the peak number (starting from one), then the approx- 
imate locations of the TT peaks are (n — 1 /4)£ s , while the EE 
peaks are at (n - 3 /4)£ s . To investigate how well our EE band- 
powers constrain the peak spacing, p hase and amplitude w e 
carry out an analysis similar to that of Readh ead et al.l ((2004). 
Subtracting the smoothed version of the LCDM EE spectrum 
from the un-smoothed (i.e. subtracting the green from the red 
curves in Figure l20i i results in a series of approximately si- 
nusoidal modulations whose fractional amplitude dies away 
close to linearly from the 2nd to the 9th peaks. This fact al- 
lows us to generate "toy-model" EE spectra which follow the 
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FIG. 19. — Bandpower deviations comparing the combined spectra to left: 
LCDM, and right: the null model. The horizontal axis is multipole number. 
The respective mo del e xpectation values have been subtracted from the spec- 
tra shown in Figure fTTI and the result divided by the error bars. Noted on each 
panel is the x 2 /ndf and the probability to exceed this value by chance (x 2 
being calculated using the bandpower covariance matrix). Note the differing 
y-axis ranges in the right column. Note also that the BB LCDM expectation 
values contain a significant leakage contribut ion an d this is why the left and 
right BB panels differ so much — see Section flO ll for details. 



envelope of LCDM using 



t(£) = S (i) + a v(£) sin ( 2tt- + p 



(21) 



where s(£) is the smoothed version of the LCDM spectrum, 
v(£) is the linear falloff, a is a re-scaling of the amplitude, and 
p is the phase. We then fit this model to the data — Figure |2T1 
shows the results. We find i s = 306 ± 10, p = 13° ± 33° and 
a = 0.86 ±0.17, consistent with LCDM, as is shown in the 
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FIG. 20. — Comparing the combined EE spectrum (blue points) to LCDM 
(red curve) and a smooth curve without peaks (green curve). The x 2 /ndf, 
and the probability to exceed this value by chance, is noted for each model 
(y 2 being calculated using the bandpower covariance matrix). 

right p art of the figure. Using WMAP TT data iPage et al.l 
(2003) found the spacing between the first and second peaks 
to be 315 ± 2. The consistency of peak phases and spacings 
between temperature data and these new QUaD EE results 
constitutes yet another confirmation of the acoustic o scilla- 
tion paradigm of CMB anisotropics. Readh ead et al.l (2004) 
allowed only the phase and amplitude to be free parame- 
ters finding p = 24 ±33° — making this restriction we find 
p = 9±13°. 
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500 1000 1500 2000 
multipole 



2000 - 

1800 

1600 ■ 

1400 ■ 

a 1200 
o 

E 1000 

3 

E 800 ■ 
600 ■ 
400 
200 ■ 
- 



LCDM TT peaks 
LCDM EE peaks 
- QUaD EE best fit 
WMAP TT peaks 




2 4 6 

peak number 



FIG. 2 1 . — Left: Fitting a "toy-model" to the combined EE spectrum to de- 
termine the peak spacing, phase and amplitude, and the uncertainties thereon. 
The red line is the initial LCDM model, which is then smoothed and sinu- 
soidally modulated with a range of peak spacings from 280 to 320 to generate 
the example family of curves shown in green. The best fitting model curve is 
shown in magenta, and the associated parameter constraints indicated. Right: 
The location of the TT and EE peaks under LCDM as compared to the best 
fit mo del and its uncertainty range. The WMAP points are from IPage et all 
<2003l) . 



9.4. BB limits and comparison to other experiments 

The BB results shown in Figure[l7]are consistent with zero 
sky power, and we therefore interpret these results as upper 
limits. To convert the observed values into confidence lim- 
its we find the 95% integral point of the positive part of the 
bandpower probability density function (which is assumed to 
be Gaussian with mean and spread as in FigurefTTb. Figurel22l 
shows a comparison of these limits, and our other spectra, to 
published results from other experiments. In this figure EE is 
shown on a log-scale, and hence to avoid clutter only band- 
powers whose center value is more than twice the distance 



between the center value and the lower end of the 68% confi- 
dence limit are shown. For TE, EE and BB QUaD breaks 
new ground — for TT the high I precision is comparable 
to ACBAR, although the beam uncertainty is larger in the 
present analysis. 

9.5. E and B maps 

In Figure[l7]it is clear that we detect dramatically more E- 
mode power on the sky than B-mode. Another way to visual- 
ize this is in the image plane. Having converted the apodized 
Q and U maps to E and B Fourier modes as described in Sec- 
tion [6] we can take the inverse Fourier transform to produce 
E and B maps. To enhance the signal to noise we apply a 
Fourier space filter 3 . Figure|23]shows the result — the E map 
contains far more structure than the B map. To confirm that 
the residual B signal is consistent with noise the figure also 
shows equivalent deck jackknife maps. 

10. SYSTEMATIC ISSUES 

There are a wide range of systematic effects which poten- 
tially mix T into pol. and/or EE into BB. In practice, as seen 
in Section |9T2l our results are consistent with the LCDM pre- 
diction of lensing BB, which is an order of magnitude smaller 
than our noise induced bandpower uncertainty, and hence ef- 
fectively zero for the purposes of our experiment. 

While there are many ways to produce false BB signal, it 
is virtually impossible that true BB power could somehow be 
canceled out through systematic effects. The fact that our BB 
results are consistent with zero is therefore powerful evidence 
that systematic mixing effects have been controlled to the re- 
quired level of precision. 

Note that from an instrumental point of view this is a for- 
tunate accident of cosmology — the theoretical prediction is 
that the CMB sky presents us with a high purity "test pattern" 
against which we can validate our experiment. If we did see 
BB power at a level greater than the LCDM prediction exten- 
sive investigation of possible systematics would be required 
to attempt to verify that it was real. But since we do not see 
any BB further investigation of mixing systematics is arguably 
unnecessary. However, for completeness, we present some 
additional discussion below. 

10.1. EE to BB mixing due to the sky cut 

Since the transform from Q,U to E,B is non-local, for 
less than full sky coverage leakage occurs from the EE spec- 
trum to BB (and vice versa, although this is irrelevant under 
LCDM). In this analysis we deal with this effec t through the 
cross-spectra BPWF's discussed in Section [6~6l For LCDM 
our BB expectation values peak at 0.7 /iK 2 for I = 360 and 
fall rapidly to higher multipoles. At the current level of sen- 
sitivity including the cross spectral BPWF's is just starting to 
become necessary — as we see from the increased y 2 value 
going from the LCDM to null model BB panels in Figure [191 

Our signal only simulations include several other effects 
which potentially produce false BB signal (including beam 
centroid mismatch and grid angle uncertainty — see Sec- 
tion 15. U . However for the standard simulation parameters 
we find that the sky cut effect is dominant since the mean BB 
simulated spectrum follows the expectation values calculated 
using the cross-spectral BPWF's. 

3 Here we use a Weiner filter assuming the LCDM EE spectrum, but a 
band pass filter 200 < I < 2000 produces a qualitatively similar result. 
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ratio of the main lobe gains which we wish to know. Sky 
dips extending over a much larger zenith angle range (5 to 
45°) follow the expected sec(6>) dependence very closely (< 
1% rms residual), indicating that the elevation nods are not 
distorted by sidelobe response. 

The two detectors of each pair share a common feed horn 
and filters, but one might worry that the bandpasses of the 
fore and aft detectors might still differ. Since the atmospheric 
emission has a different frequency spectrum from the CMB 
this might lead to a systematic error in the relative gain within 
each pair. 

To test how such errors would play out in practice we ran 
some special simulations where the input sky maps contain T 
only, and where the detector pair gains are deliberately sys- 
tematically distorted such that gf ore /gaft = 1-03 — far worse 
than the < 1 % constraint on possible mismatch which we de- 
rive from RCW38 observations. We find that the resulting EE 
and BB spectra follow the same form as the input TT spec- 
trum with a peak of 0.4 /iK 2 at I = 200. It is important to note 
that even for a systematic error like this there is still consider- 
able averaging down as pairs of different angles, at different 
telescope orientations, leak T into, for instance, both +Q and 
—Q within a given map pixel. Making the gain ratio errors a 
random 3% across the focal plane (but fixed over time) the 
averaging down is much more effective and the peak leakage 
becomes 0.1 /iK 2 . 



FIG. 23. — E and B signal and (deck) jackknife maps. The color scale is 
±20 /iK. These maps have been apodized and filtered to enhance signal-to- 
noise — see text for details. 

The original MASTER approach dHivon et all 120021) in- 
cludes a "de-mixing" correct ion, and this was subsequently 
extended to polarization (e.g. lBrown et al.ld2005l) ). Such tech- 
niques reduce spectral mixing in the mean, and hence also re- 
duce the cross-sp ectra BPWF's. More recen tly a technique 
was proposed by ISmith & Zaldarriagal (120071) which results 
in much lower mixing from E to B within any given realiza- 
tion — we have implemented this in the flat sky case, and 
confirmed with simulations that it works. However since this 
complication is unnecessary we do not include it in this paper. 

10.2. Curved sky versus flat sky 

The analysis presented here uses flat sky power spectrum 
estimation while most recent CMB analyses have been con- 
ducted in the spherical harmonic basis. By using full curved 
sky maps as the input to our signal timestream simulation, our 
power estimation is normalized to recover the curved-sky in- 
put power. In addition these simulations empirically test for 
problems associated with these flat-sky estimators, which are 
found to be negligible for our small (~6x 6°) patch of sky. 
As mentioned above EE to BB leakage is dominated by the 
cut sky effect and is well reproduced in the simulations by the 
cross-spectral BPWF's — for QUAD the use of flat-sky power 
estimators does not contribute significantly to this leakage. 

10.3. Relative Gain Calibration 

In Section 13.21 the elevation nod based method we use to 
equalize the detector gains was described. This method ap- 
pears to be very accurate, and we do not include any scatter 
in the standard simulations. However, since error in the pair 
relative gains leads to T to pol. leakage there are a couple of 
issues which one might worry about. 

The elevation nods integrate over the full beam including 
any far sidelobes, whereas when mapping the CMB, it is the 



10.4. Pair beam mismatch 

As mentioned in Section 15.11 we observe repeatable beam 
centroid offsets between the two halves of each detector pair, 
fixed in the instrument frame, with rms magnitude of ~ 0.1'. 
T to pol. leakage introduced by this effect is included in 
the standard simulations but makes a negligible contribution. 
However for interest we run some special simulations with 
only T input introducing random pair centroid offsets with 
rms magnitude of 2' — twenty times the observed value. We 
find that the resulting leakage has a broad peak at ~ 2 /iK 2 
around I ~ 1000 for 100 GHz, and - 1 /iK 2 for 150 GHz. 
In addition we run simulations under the totally unrealistic 
scenario that all pairs are systematically offset in the same di- 
rection by 1' and find leakage of 6 /iK 2 at £ ~ 1000. 

For the standard simulations we use the individual chan- 
nel major/minor fit widths and orientation angles as measured 
using PKS0537-441. Some of the apparent variation between 
channels in these observations is measurement noise, but there 
is some real variation, and < 10% elongation. To investigate 
the impact of differential beam size we run special T only sim- 
ulations with both the major and minor axis FWHM for one 
half of each pair systematically inflated by 1 .4' — the result- 
ing leakage peaks at ~ 1 .5 /iK 2 above £ ~ 1500. For differen- 
tial elongation, where we inflate only one axis of one half of 
each pair (with a common orientation angle for all pairs), we 
find - 1 /iK 2 above £ ~ 1500. 

10.5. Polarization angle 

As seen in Section 19.21 our results are consistent with 
LCDM — we detect considerable EE power and no BB. This 
is in itself confirmation that the assumed polarization angles 
of the detectors are known to sufficient accuracy. If we did 
see significant BB we might suspect that it was false signal 
due to incorrect angles. But it would be nearly impossible for 
a sky which truly had B-mode power to appear not to due to 
incorrect detector angles. 
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To confirm this we re-generated the real maps using detec- 
tor angles systematically biased from the best estimate va lues 
by far more than the 1 ° estimated uncertainty (see Section l4~3l 
and the Instrument Paper). For a 5° bias there is almost no 
effect — the x 2 PTE versus LCDM for the BB spectrum (as 
shown in Figure [191) falls from 0.56 to 0.21. Only at 10° do 
we start to see significant extra signal in the BB bandpowers, 
and x 2 failure versus LCDM for BB (and EB). 

10.6. Polarization efficiency 

As mentioned in Section 14.31 and our Instrument Paper, 
we measure the polarization efficiency of our detectors to be 
e = 0.08 with rms scatter of 0.015. An additional calibration 
factor 7 = (1 — e)/(l +e) is then applied to the pair difference 
data. Random errors in e will average down, while a system- 
atic error will translate into a shift in the absolute calibration 
of the polarization power spectra ~ 4 times as large (includ- 
ing the additional doubling when going from units of temper- 
ature to power). Since we estimate the uncertainty on e to 
be < 0.02 the implied uncertainty on the polarization spec- 
tra is sub-dominant to the overall absolute calibration uncer- 
tainty. 

10.7. Moon pickup 

QUaD has a variety of far sidelobes as described in the ac- 
companying Instrument Paper. It is not clear which of these 
produces the bulk of the highly polarized ground pickup men- 
tioned in Section [331 — however since the field differencing 
is so effective at removing this contamination this is probably 
of academic interest only. However, any source of contam- 
ination which moves with respect to the ground will not be 
removed by field differencing. 

Part of our basic low level data reduction infrastructure in- 
volves making single pair sum and difference maps for each 8 
hour block of observations. On certain days when the Moon 
is high above the horizon, but at a very large angle from the 
telescope pointing direction, we see obvious stripes in these 
maps. Even after cutting these visibly contaminated days 
from the analysis, we saw a strong peak towards zero in the 
PTE distribution as show in Figure [14] and some x 2 values 
were much too large. 

Extensive effort has been required to elucidate the coupling 
mechanism by which the Moon enters the CMB field data. 
We have determined that radiation reflects off the inside sur- 
face of the foam cone creating a polarized ring sidelobe at 
sw 100° from the main beam. The shape of the pickup across 
any given scan has been successfully modeled using the posi- 
tion of the Moon relative to the telescope pointing direction, 
the feed offset angle, and the detector polarization angle (see 
the Instrument Paper for details). 

Using this model we cut data periods which are potentially 
contaminated by the Moon. For this analysis the cut is simple 
and quite aggressive — for any day where the Moon passes 



within a generous band about the ring sidelobe we simply 
discard the entire day. A future analysis could retain slightly 
more data by selectively cutting channels, and using a time 
granularity shorter than a whole day. Performing this cut we 
reject an additional 59 days of observation but as seen in Sec- 
tion|7]the jackknife tests pass, meaning that we can be confi- 
dent of the final power spectrum results. 

11. CONCLUSIONS 

We have described the observations, data reduction, simu- 
lation and power spectrum analysis of the QUaD experiment. 
The results reported here are from 143 days of data taken in 
the second and third (final) seasons of observation, employing 
a conservative Moon cut, and simple lead-trail differencing. 
A future analysis may be able to include more data, and/or 
reduce the information lost in the ground removal. 

The three sets of power spectra, 100 GHz, 150 GHz and 
frequency-cross, have been subjected to an extensive set of 
jackknife tests and residual systematic contamination has 
been shown to be undetectable above the instrumental noise. 

The combined spectra improve very considerably in sensi- 
tivity over previous results, and are consistent with LCDM — 
the standard cosmological model has passed yet another strin- 
gent test. Furthermore we find that a smooth curve is a very 
poor fit to the observed EE spectrum — acoustic peaks in the 
EE spectrum have been detected with high significance for 
the first time. The impact of possible instrumental systematics 
has been considered in detail, but in fact the tight upper limits 
on BB power obtained are in themselves extremely powerful 
evidence that such effects are adequately controlled. 
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